The transition temperature of the dilute interacting Bose gas 
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We show that the critical temperature of a uniform dilute Bose gas increases linearly with the s-wave 
scattering length describing the repulsion between the particles. Because of infrared divergences, 
the magnitude of the shift cannot be obtained from perturbation theory, even in the weak coupling 
regime; rather, it is proportional to the size of the critical region in momentum space. By means of 
a self-consistent calculation of the quasiparticle spectrum at low momenta at the transition, we find 
an estimate of the effect in reasonable agreement with numerical simulations. 



Determination of the effect of repulsive interactions on 
the transition temperature of a homogeneous dilute Bose 
gas at fixed density has had a long and controversial his- 
tory [1-5]. While |IJ predicted that the first change in 
the transition temperature, T c , is of order the scatter- 
ing length a for the interaction between the particles, 
neither the sign of the effect nor its dependence on a 
has been obvious. Recent renormalization group studies 
j| predict an increase of the critical temperature. Nu- 
merical calculations by Griiter, Ceperley, and Laloe ||, 
and more recently by Holzmann and Krauth 0, of the 
effect of interactions on the Bosc-Einstein condensation 
transition in a uniform gas of hard sphere bosons, and 
approximate analytic calculations by Holzmann, Griiter, 
and Laloe of the dilute limit ||, have shown that the 
transition temperature, T c , initially rises linearly with a. 
The effect arises physically from the change in the energy 
of low momentum particles near T c [|[ . Here we analyze 
the leading order behavior of diagrammatic perturbation 
theory, and argue that T c increases linearly with a. We 
then construct an approximate self-consistent solution of 
the single particle spectrum at T c which demonstrates 
the change in the low momentum spectrum, and which 
enables us to calculate the change in T c quantitatively. 

We consider a uniform system of identical bosons of 
mass m, at temperature T close to T c and use finite tem- 
perature quantum many-body perturbation theory. We 
assume that the range of the two-body potential is small 
compared to the interparticle distance n -1 / 3 , so that the 
potential can be taken to act locally and be character- 
ized entirely by the s-wave scattering length a. Thus we 

1/2 

work in the limit a <C A, where A = (2irh /rnksT^ 
is the thermal wavelength. (We generally use units 
h = k B = 1.) 

To compute the effects of the interactions on T c , we 
write the density n as a sum over Matsubara frequencies 
u) v = 2-nivT {y = 0, ±1,±2, ...) of the single particle 
Green's function, G(k,z): 
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where 



k 2 

G^ 1 (k, z) = z + n S(fc, z), 

2m 
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with /x the chemical potential. The Bose-Einstein con- 
densation transition is determined by the point where 
G- x (0,0) = 0, i.e., where £(0,0) = /i. 

The first effect of interactions on £ is a mean field 
term £ m / = 2gn, where g = 4Trh 2 a/m; the factor of 
two comes from including the exchange term. Such a 
contribution, independent of k and z has no effect on the 
transition temperature, as it can be simply absorbed in a 
redefinition of the chemical potential. To avoid carrying 
along such trivial contributions we define: 



= -(fj>- 2gn)- 



(3) 



The quantity £ may be interpreted as the mean field cor- 
relation length. In the mean field approximation, £ be- 
comes infinite at T c ; however, in general, it remains finite, 
and functions here as an infrared cutoff. 

Because the effects of interactions are weak, one could 
imagine calculating the change in T c in perturbation the- 
ory. However, such calculations are plagued by infrared 
divergences. Power counting arguments reveal that the 
leading contribution to the self-energy, £(fc <C C _1 ; 0), of 
a diagram of order a n has the form: 



£n 



(4) 



In perturbation theory about the mean field, with the 
mean field criterion for the phase transition, £ — > 00 at 
T c , all £„ diverge, starting with a logarithmic divergence 
at n = 2. More generally, the approach of C towards 
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A 2 /a in magnitude signals, according to the Ginsburg 
criterion, the onset of the critical region. Beyond, per- 
turbation theory breaks down, since all S n in Eq. (Q) are 
of the same order of magnitude. 

Even though the theory is infrared divergent, we can 
isolate the leading correction to the change in T c , which, 
as we show, is of order a. Since the infrared divergences 
occur only in terms with zero Matsubara frequencies we 
separate, in Eq. (|l|), the contribution of the v = terms, 
writing 



/d 3 k 
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where G iy ^o(k) is the sum of terms with v ^ 0. Similarly 
the density of a non-interacting system with condensa- 
tion temperature T is given by 



n°(T) 



-T 
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(G° = o(fc) + G°^ (fc)) 
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where C(3/2) = 2.612. Since the non-zero Matsubara 
frequencies regularize the infrared behavior of the mo- 
mentum integrals, the dependence of the v ^ terms in 
Eq. (|J) on a is non-singular at T c . These terms, of order 
a 2 at least, can be neglected. Thus to order a, 



n(a,T c )-n°(T c ) 



d 3 k 



(G (k)-G° (k)). (7) 



To calculate the change in T c at fixed density we 
equate n(a,T c ) at T c with n°(T c °) at the free particle 
transition temperature T° and observe that n°(T c ) = 
(T c /T°) 3 / 2 n°(T°); thus in lowest order the change in 
transition temperature AT C = T c — T c ° is given by 



n°(T c °) = T c j (GWo(fc) - G u=0 (k)). (8) 



3AT C n/ nN f d 3 k 

where AT r = T, - T° Thus 
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U(k) 



37rC(3/2)7 k 2 + U{kY 
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where U(k) = 2m(S ly=0 (fc) — /i). 

Equation (^) for the leading correction to the critical 
temperature is crucial. The criterion for spatially uni- 
form condensation is that U(0) = 0; above the transition, 
U(0) > 0. At the transition, k 2 + U{k) > for k > 0. At 
large wavenumbers, U — > l/( 2 > 0, and in the critical re- 
gion, as we discuss below, U is also positive. Although we 
have not proved it rigorously, numerical simulations indi- 
cate that U is generally positive for k > 0, which implies 
that the integral in Eq. ([)]) and hence AT C is positive. 

In the critical region, k < k c , where k c defines the 
scale of the critical region in momentum space, G„ = o 



has the scaling form § G~i {k) = -k 2 ^ ] k?F '(fc£); £ is 
the coherence length which diverges at T c as |T — T c \~ v , 
and F is a dimensionless function, with F{po) ~ 1. The 
critical index, 77, is given to leading order in the e = 4 — d 

expansion by e 2 /54 Q. At T c , G~l (k) k 2 ^k^, 

and thus U ~ +fc 2_r; . Both terms in Eq. (Q) give a 
contribution of order k c , so that AT C /T C ~ fc c . As we 
shall see, fc c ~ a/ A 2 , and hence AT C /T C ~ a/ A. 

To study the leading behavior in a quantitatively, we 
need concentrate only on the v = sector where the 
full finite temperature theory reduces to a classical field 
theory pot] defined by the action: 



S{4>(r)} = 



2mT 



( V0*(r) • V0(r) 



■^(r)] 2 +4M\Hr)\ 2 ~ Mr)\ 2 )) : 



(10) 



the probability of a given field configuration entering 
the computation of expectation values, is proportional 

The classical theory is ultraviolet divergent, but su- 

perrenormalizable. The divergences appear only in the 

(2) 

two-loop self-energy, S^ - effectively the second order 
self-energy written in terms of the full G„ = o rather than 
the zeroth order Green's functions - and can be removed 
by simple renormalization of the mean field coherence 
length, Since, henceforth, we calculate only in the 
classical theory, we drop the subscript v = . The sec- 
ond order self-energy is 



(id 



where eu — (k 2 + C 2 )/2m, and the {y = 0) particle-hole 
bubble, 



B(q) 
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is given by 



B(q) 



T\ 4 



(12) 
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b(x) -> 1/x for x > 1 and 6(0) = 

The integral in Eq. ([ll]) is logarithmically divergent 
in the ultraviolet. But in the full theory the momen- 
tum integrals are cut off by distribution functions, / — 
(e k / 2mT — and the ultraviolet behavior is regular. 

To control this divergence we introduce an ultraviolet 
momentum cutoff, A, in the classical theory, recognizing 
that it is in fact effectively determined in the full theory. 
Then 



2mS(fc) 



a 2 f H 

327T 2 — / dxxb(x)L(k(^,x), 
A Jo 



(14) 
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where 



(15) 



The divergent part of the integral comes from the large x 
tail of b(x), and contributes -128(a/A 2 ) 2 ln(AC) to 2mE. 

More generally we carry out a diagrammatic ex- 
pansion of £ in terms of the self-consistent v = 
Green's function, defined by 2mG~ 1 (k) = —k 2 + £~ 2 + 
2mE(fc, a, G, A). Note that the dependence of E on f en- 
ters only through the dependence of £ on G. We define 
a renormalized mean field coherence length by 



1 

T 2 

^R 



1 



128 



In(ACfl) 



Then G 1 (k) is given by 



2mG~ 1 (k) = k 2 + Q R 2 + 2m£ F (fc, a, G), 



(16) 



(17) 



where 



£ F (fc, a, G) = £(fc, a, G, A) + 128 



A 5 



ln(AC fl ) (18) 



is independent of A. As a function of £ F , the Green's 
function is independent of the cutoff. 

In fact, a simple power counting argument shows that 
the finite part of the self-energy has the form 



£ F (fc, a, G) 



1 



2mC 



2 &(Kr, J) 



where 



J = <r/\ 2 . 



(19) 



(20) 



To see this structure we note that a term in the self- 
energy of order a n is the product of a dimensionless func- 
tion of k(,R times the E n of Eq. ([|) , with ( replaced by 

Cr & " 

The criterion for condensation, ( R 2 + 2m£ F (0, a, G) = 
0, implies that 



1 + <r(0, J) = 0. 



(21) 



Since cr(0) is a well-behaved function of only the param- 
eter J, Eq. (|l] determines the critical value of J = J* 
for condensation, a dimensionless number independent 
of the parameters of the original problem. At conden- 
sation, the renormalized mean field coherence length £r 
tends to infinity as a — > 0, with the product a^R fixed, 
thus preventing a perturbative expansion in a. 

At condensation U(k) = (cr(k( R , J*) + 1)/£ R) and 
Eq. (g) implies the change in T c 



T, 



1 



a ( 4 
A V3ttC(3/2) J* 



dx- 



cr(x, J*) + 1 
' x 2 + cr(x, J*) + 1 



(22) 



Since J* is determined by the condition (pl|), the result 
for AT C /T C is linear and expected to be positive in a/ A. 

We turn now to calculating AT C explicitly within a sim- 
ple self-consistent model based on taking only the zero 
frequency component of the leading two-loop approxima- 
tion self-energy, given by Eq. (|ll]). We construct the e p 
as self-consistent quasiparticle energies at the transition, 
i.e., solutions of the equation: 



G-\k,e k ) = = e k 



2 m 



(S(fe)-S(O)). 



(23) 



The low momentum behavior of e k is determined by 
a familiar argument Jl2]| . In order that the integral ( |Tl| ) 
converge in the infrared limit, e k must behave, modulo 
possible logarithmic corrections, as ~ k a , where a < 2. 
In this case, the term k 2 /2m in fl23| ) can be neglected 
at small k. We then expand about k = 0. For 

1 < a < 4/3 the self-energy is sufficiently convergent 
that S(fc) — E(0) ~ fc 2 at small k, and thus cannot be 
the correct self-consistent solution. For a with 4/3 < 
a < 2 one has £(fc) - S(0) - +fc 6 ~ 3a , so we find self- 
consistency, E(fc) - E(0) ~ fc a , for a = 3/2. We write 
the small k part of the spectrum as 



efc 



fc c 1 / 2 fc 3 / 2 /2m. 



(24) 



Here k c is the wavevector around which the k z l 2 at low 
k crosses over to the k 2 /2m free-particle behavior. 

To extract the low momentum structure, below the 
scale fc c , we evaluate the most divergent terms of 



SW - S( o, = -Vr/|L B („(-l--I 



(25) 



at small k. Since the the A; 3 / 2 structure arises from the 
small q behavior of the integral; we evaluate the bubble 
B(q), Eq. (O), at small q with the spectrum ( p4[ ) for 
k < k c and Ar/2m for k > k c . Then 



5(9) 



4m 
7rA 2 fc, 



■(ln(fc c /g)+c), 



with c « 2 + 2 In 2 - tt/2 =1.816. Thus, 
1024tt / a \ 2 / k 



3/2 



(26) 



(27) 



Identifying the right side of Eq. @ with kl /2 k 3/2 /2m, 
we derive 



(28) 



As expected, the scale of the unusual low momentum 
structure is a/ A 2 . 

Let us, for a first quantitative estimate, assume a 
spectrum at T c of the form e k = kl^ 2 k 3 ^ 2 /2m for 
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fc <C k c , and (fc 2 + fc 2 )/2m for k 3> k c . We smoothly 
interpolate between these limits, writing U (fc) = 
fc c 1/2 fc 3 / 2 / (1 + (fc/fc c ) 3 / 2 ). Thus JdkU/(k 2 + U) ~ 1.2fc c , 
so that with Eq. Q, 



AT C 



2.9an 1/3 . 



(29) 



By comparison, Griiter, Ceperley, and Laloe Q find 
AT C /T C rj 0.34a?! 1 / 3 , while the more recent calculation of 
Holzmann and Krauth yields AT C /T C « (2.3±0.25)an 1 / 3 . 
The agreement of the numerical coefficient, given the 
simplicity of the approximations in evaluating the effect 
of interactions on the transition temperature, is satisfy- 
ing. As will be reported in a fuller paper Q, this es- 
timate agrees with that derived from the numerical self- 
consistent solution of Eq. (|2^) . 

The lowest two-loop calculation does not account fully 
for the modification of the transition temperature; in- 
deed, at the critical point, all diagrams become compa- 
rable |l^,|l^|. Consider, for example, summing the bub- 
bles describing the repeated scattering of the particle- 
hole pair in B Q , thus replacing B in Eq. ([ll]) by 



B cS {q) = 



B{q) 



l + 2gB(qY 



(30) 



where the two accounts for the exchange terms. The 
denominator at small q, from Eq. (Eq), is given by 



1 + 2gB(q) = 1 



32a 

Vkr. 



(ln(fc c /q) + c). 



(31) 



Since k c ~ a/ A 2 , the correction is of order unity, and 
serves to modify the spectrum, recalculated from Eq. ( |25| ) 
with (|H), from fc 3 / 2 to fc 2 ^, with @ r? ~ (l/2)-l/(2c+ 
fc c A 2 /16a) ~ 0.36. 

To estimate J*, we calculate S(0) from Eq. ([ll]) with 
the 3/2 spectrum and the leading log in B(q), Eq. (poj), 
and neglect the contribution for q > k c . Then E(0) ~ 
-K 2 a 2 /2mA 4 , and cr(0, J) ~ —k 2 J 2 , so that at T c , J* ~ 
1/k = 3/(32-^(2 + 3c)). The self-consistent solution of 
Eq. (||) yields Q J* ~ 0.07. 

The modification at T c of the spectrum of particles 
at low momenta should have direct experimental conse- 
quences in trapped Bose condensates. While a fc 2 /2m 
particle spectrum yields a flat distribution v 2 dn/dv of 
velocities, a more rapidly rising spectrum, e.g., the fc 3 / 2 
discussed here, depletes the number of low momentum 
particles. These effects become more pronounced with a 
larger number of particles and flatter traps, as the level 
spacing ceases to control the low energy behavior. 
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